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In this work we derive a general formula for the charge pumped in a superconducting nanocircuit. 
Our expression generalizes previous results in several ways, it is applicable both in the adiabatic and 
in the non-adiabatic regimes and it takes into account also the effect of an external environment. 
More specifically, by applying Floquet theory to Cooper pair pumping, we show that under a 
cyclic evolution the total charge transferred through the circuit is proportional to the derivative of 
the associated Floquet quasi-energy with respect to the superconducting phase difference. In the 
presence of an external environment the expression for the transferred charge acquires a transparent 
form in the Floquet representation. It is given by the weighted sum of the charge transferred in 
each Floquet state, the weights being the diagonal components of the stationary density matrix of 
the system expressed in the Floquet basis. In order to test the power of this formulation we apply it 
to the study of pumping in a Cooper pair sluice. We reproduce the known results in the adiabatic 
regime and we show new data in the non-adiabatic case. 

PACS numbers: 85.25.Cp, 03.65.Yz, 03.65.Vf 



I. INTRODUCTION 

In a mesoscopic conductor a dc charge current can be 
obtained, in the absence of applied voltages, by cycling in 
time two (or more) external parameters, e.g. gate volt- 
ages and/or magnetic fluxes, which govern the transport 
properties of the systenffl Adiabatic charge pumping 
refers to the regime when the variation of the external 
parameters is slow as compared to the characteristic time 
scale of the system. 

In the scattering approach to quantum transport the 
pumped charge in an adiabatic cycle can be expressed 
in terms of derivatives of the scattering amplitudes with 
respect to the pumping parameters^. In the opposite 
regime of Coulomb blockade, with several metallic islands 
connected to each other by small tunnel junctions, a pe- 
riodic modulation of the externally applied gate- voltages 
leads to a periodic lifting of Coulomb blockade, thus en- 
abling the transfer of exactly one electron per period 
through the device. Experimental evidence for para- 
metric charge pumping in normal metallic systems in 
the regime of Coulomb blockade has been obtained in 
Refs. 13141 Over the last decades charge pumping has at- 
tracted the interest of many research groups working on 
very different aspects of this phenomenon ranging from 
its metrological applications to its intimate relation with 
the fundaments of quantum theory (see Ref. [5] and refer- 
ences therein). 

Originally motivated by the aim of achieving quantized 
charge pumping in the GHz range, a great deal of atten- 
tion has been devoted in the last two decades to super- 
conducting systems. The first experiment in this context, 
performed by Geerligs et aP, showed that the degree of 
quantization of the pumped charge was not as good as in 
the normal case. As it was later discussed by Pekola et al 
7 , the main source of inaccuracies is related to the overall 
coherence of the superconducting system. This apparent 
disadvantage (appropriate designs of the superconduct- 



ing circuit may overcome this difficulty) however turned 
out to be a precious source for the investigation of fun- 
damental properties of quantum theory in macroscopic 
systems. 

If only superconducting leads are present, at low 
enough temperature, pumping is due to the adiabatic 
transport of Cooper pairs. Besides the dependence of 
the pumped charge on the details of the cycle, in su- 
perconducting pumps there is an additional dependence 
on the superconducting phase difference since the overall 
process is coherent. Cooper pair pumping has been thor- 
oughly investigated^^ in the last d ecades. In a series of 
experiments the Helsinki g r0U pU2im has shown the coher- 
ent properties of Cooper pair pumping and, very impor- 
tantly, provided the first experimental demonstration of 
the relation between Cooper pair pumping and the Berry 
phase acquired by the system during its cyclic evolution. 
A connection between Berry phase and pumped charge 
in superconducting nano-circuits has been already estab- 
lished theoretically in Refs. I9I12I13I (see also Refs. 1191201 
where this relation was found for mesoscopic normal con- 
ductors). 

Berry phases in macroscopic systems such as supercon- 
ducting circuits have been studied in Refs. I13l21lj2"3l and 
very recently experimentally demonstrated by the ETH 
groupJ^. The large body of theoretical understanding 
and the spectacular experimental control which lead to 
unveil the coherent properties of pumping and its rela- 
tion to geometric phases are important steps toward s the 
implementation of geometric quantum computation^^!] 
with superconducting devices. 

While a lot has been found in the relation between 
geometric phases and pumping in closed quantum sys- 
tems, the role of an external environment constitutes, 
with the notable exception in few very recent papers 27 29 , 
an almost unexplored territory. The study of geomet- 
ric phases in the presence of decoherence and dissipation 
has started only recently with few exceptions though, cer- 
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tainly prompted by the interest in quantum computation 
(see for example Ref . l30|) . Together with many features 
common in the theory of open quantum systems, the 
analysis of decoherence in geometric interfcrometry rises 
several distinct issues that are of interest both as funda- 
mental questions in quantum mechanics and in quantum 
computation. The adiabatic evolution, for example, can- 
not occur arbitrarily slow, as decoherence would destroy 
any interference. This implies that the decoherence pro- 
cesses should be analyzed in close connection with non- 
adiabatic corrections^. 

Being Cooper pair pumping a geometric quantum ef- 
fect, it is natural to ask oneself to which extent an ex- 
ternal environment modifies its characteristics. Not only, 
this question is relevant for a detailed comparison with 
experimental data where an external bath is unavoidably 
present, but also it may shed additional light in the role 
of dissipation on geometric quantum phenomena. It is 
not a priori obvious, for example, that a relation between 
pumping and Berry phases (of any sort), will survive in 
an open system. As already mentioned, till now this 
problem was tackled in Rcfs. l2"TH2l?l where a generalized 
master equation to consistently account for the combined 
action of the driving and dissipation was derived. Appli- 
cation to the Cooper pair sluice^ showed that in the zero 
temperature limit the ground state dynamics, and con- 
sequently pumping, is not affected by the environment. 

Stimulated by the results obtained in Refs. |2"7H2T)1 in 
this paper we further investigate the relation between 
pumping and geometric phases. We derive an expression 
for the pumped charge which generalizes previous results 
in several aspects. It is valid also under non-adiabatic 
conditions and in the presence of an external environ- 
ment. The key to our approach is to apply Floquet the- 
ory to Cooper pair pumping. We will show that under 
cyclic evolution of the system the total charge transferred 
through the circuit is proportional to the derivative of the 
associated Floquet quasi-energy with respect to the su- 
perconducting phase difference (a result which is valid 
also in the case of a non-adiabatic evolution). In the 
presence of an external environment the expression for 
the transferred charge is easily generalized in the Flo- 
quet representation. It is given by the weighted sum of 
the charge transferred in each Floquet state, the weights 
being the diagonal components of the stationary density 
matrix of the system expressed in the Floquet basis. The 



central result of our work is Eq. ( 20 ) ; it embraces all the 



limits considered so far in the literature and allows to 
investigate new regimes. Furthermore it suggests the use 
of a series of well known numerical schemes to compute 
the pumped charge. 

The paper is organized as follows. In the Sections |TT| 
and [HI] we will introduce the basic ingredients needed in 
the derivation of the pumping formula. In Section [IT] we 
formulate the problem of Cooper pair pumping in super- 
conducting circuits while in Section |III| we provide the 
necessary tools of the Floquet theory of driven quantum 
systems both in the closed and open cases. The formula 



for the pumped charge will be derived in Section IV Here 
we will discuss in which aspects our results generalizes 
previous works. As an example we will apply our deriva- 
tion to the Cooper pair sluice which was experimentally 
realized in the Helsinki group. In Sectio n [V] various dif- 
ferent limits will be discussed. Section IVII will contain 
the conclusions of the present work. 



II. COOPER PAIR PUMPING 

Charge pumping in Josephson networks consists in a 
coherent periodic manipulation of the collective state of 
the Cooper pairs in the array. In this Section wc define 
the setting and review, for later convenience, the relation 
between Cooper pair pumping and Berry phases. 

A Cooper pair pump consists of a Josephson network 
connected through Josephson junctions to two supercon- 
ducting leads (see a sketch of the setup in Fig[TJ). The 
system is phase biased, i.e. the two superconducting elec- 
trodes are kept at a finite phase difference tp — tpn — tp L 
where Pr/l is the phase of the superconducting order 
parameter of the right/left lead. The Cooper pair pump 
operates by changing adiabatically in time some external 
parameters such as gate voltages, to tune the charging 
energies, or magnetic fluxes, to vary the effective Joseph- 
son couplings. We will label this set of external param- 
eters by the vector A(t) = {V g i(t),$i(t)}. In the ab- 
sence of an external environment charge transport is a 
purely coherent phenomenon. The Hamiltonian of the 
pump depends on the superconducting phases of each 
island of the network tpi (i = 1, . . . N), its conjugate mo- 
menta (i.e. the charge on each island n,-), the phase dif- 
ference across the pump and all the external parameters, 
H(t) = H Upx, (pjy; n%, n^; A(t), tp] . The state of 

the system is denoted by = X(t), tp)). 

By changing the control parameters in time a charge 
Q( tr ) will be transferred from the left to the right elec- 
trode. The total charge transferred through the pump 
during the period T is given by: 



Q 



(tr) 



f 



-2e, 



(i) 



Assuming that there are no degeneraci es in the spectrum 
and in the adiabatic limit, it was 

showrPMlthat Q(tr) 

can 

be expressed in terms of the total phase accumulated by 
the system after the cycle 



)M Q lD Q lE 



2e 



dtp dtp 



(2) 



In the previous expression jd/b are the dynami- 
cal/geometric contribution to the accumulated phase. 
The dynamical contribution is the charge transferred 
through the circuit due to the supercurrent flow. The ge- 
ometric contribution is the pumped charge. In the latter 
case the charge is an even function of the superconduct- 
ing phase difference while the contribution due to the 
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FIG. 1: A generic setup of a Cooper-pair pump. Two su- 
perconducting leads, kept at a phase bias ip = ipn — <Pl, are 
connected through a Josephson network. The system is op- 
erated in a regime where quantum effects are important. In 
order to have a pumped charge some external parameters (e.g. 
gate voltages or magnetic fluxes) are varied periodically. 



supercurrent flow is odd in ip. The different symmetry 
under reflection of the phase bias is fundamental for the 
experimental detection of the pumped charge. A gener- 
alization to the non-abelian case, i.e. in the presence of 
degeneracies in the spectrum, has been given in Ref. 1151 

In order to investigate non-adiabatic corrections and to 
generalize this result to the dissipative case, it is useful 
to re-express the pumped charge using the Floquet for- 
malism. In the next sections, after introducing the basic 
definitions of Floquet theory, we will find an expression 
for the pumped charge which in the adiabatic limit and in 
the absence of an external environment reduces to Eq. (12]) 



III. FLOQUET THEORY 

As it will become clear in the continuation of the paper, 
Floquet formalism is naturally suited to study Cooper 
pair pumping. It is important to stress already now that 
this is not a mere re-formulation of what has been done so 
far. We will show that, on the contrary, Floquet approach 
treating on the same footing adiabatic and non-adiabatic 
regimes, provides a transparent and general expression 
for the pumped charge in the case in which the super- 
conducting network is coupled to an external environ- 
ment. A Floquet scattering theory has been developed 
by Moskalets and Buttiker to study pumping in meso- 
scopic conductors^. Here we employ Floquet approach 
to study Cooper pair pumping. 

In the next sections we introduce the necessary ingredi- 
ents of Floquet theory and its use in quantum dissipative 
systems. The presentation follows closely Refs. 1331341 



A. Basics of Floquet theory 

Given a system whose dynamics is governed by a pe- 
riodic Hamiltonian H (t) — H(t + T), Floquet theorem 
states that solutions to the Schrodinger equation exist 
which have the (Floquet) form 

!*»(*)) = e-^ h \$ a (t)) , (3) 

where the state |<E> a (i)) is called Floquet mode and it is 
periodic in time (|3> Q (t + T)) = | <&„(£))) and the corre- 
sponding quasi-energy e a is real and unique up to multi- 
ples of Ml, with SI = 2ir/T. There are as many distinct 
such solutions as the dimension of the Hilbert space H. 
These solutions are linearly independent and form a ba- 
sis of the Hilbert space. An eigenvalue equation for the 
quasi-energy e a can be obtained by defining the operator 
H(t) = H(t) - ihd t 

H(t)\$ a (t)) = e a \$ a {t)) . (4) 

The Floquet modes |$ Q ,„(£)) = | exp(— inClt) with 

integer n lead to a solution identical to the one given 
in Eq.pl), but with shifted quasi-energy e a — > e ai „ = 
€ a — nn&l; hence the eigenvalues {e Q } can be mapped in 
a first Brillouin zone obeying to —Ml/2 < e < Ml/2. 

For the Hermitian operator H(t) it is convenient to 
introduce the composite Hilbert space-^ "H <X> T made by 
the tensor product of the Hilbert space H of the vectors 
representing the state of the system and the space T of 
the periodic functions in t with period T = 2tt/CI. In 
the space of the periodic functions of t we have a basis 
of Fourier vectors {exp(— mQi)} which are orthonormal 
with respect to the scalar product given by 

(m,n) = i [ T (e- mm )*e- mm dt = S m , n (5) 
1 Jo 

We define |Z) = e~ tlnt as vectors in T. We can extend 
the scalar product on % to a scalar product on % ® T 
defining 

((*i|* 2 »^^ T (*i(<)|* 2 (t))dt . (6) 

B. Floquet states and geometric phases 

The Floquet quasi-energies are intimately connected 
to geometric phases (in the present work we are inter- 
ested only in the case in which their spectrum is non- 
degenerate). Indeed they are, up to a multiplying factor, 
the phases of the eigenvalues of the evolution operator 
Us(t + T,t). The phase is defined up to 2n7T, hence an 
eigenvalue of Us (t + T, t) corresponds to infinite Floquet 
exponents obtained through translations of 2mrh/T. 

Noticing that Floquet states follow a cyclic evolution 
in the projective Hilbert space (we call C the closed path 
followed in the projective Hilbert space), it is possible 
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to express the Aharonov-Anandan geometric phase ac- 
quired during a cyclic evolution starting in the Floquet 
eigenstate a by 

1aa{C) = + l£(<£ a (t)\H(t)\<S> a (t))dt (7) 

An equivalent expression, perhaps more useful in the 
computation is 



1AA 



(C) = 27ry^ k(c ayk \c a , k ) 



(8) 



where we used the Fourier expansion of the Floquet 
modes 



IM*)> = E \ c «> 



am 



(9) 



/= — oo 



In the adiabatic limit the quasi-energy corresponding 
to the n-th eigenstates can be expressed in terms of the 
dynamic, jd, and geometric, 7c phases: 



e n = - ji[7D,n(T) +jG,n(C)] 



(10) 



C. Floquet-Born-Markov Master equation 



The Floquet basis is particularly useful to write the 
Master equation governing the dynamics of the reduced 
density matrix of a driven system when in contact with 
an external environment. We consider below the case in 
which the Born-Markov approximation is applicable. De- 
tails and subtleties of the derivation of a Master equation 
in this case are described in Ref. [331 here we merely state 
the end result for the Master equation which will be later 
used to derive a formula for the pumped charge. 

Given a quantum system interacting with an ex- 
ternal reservoir, the Hamiltonian describing the sys- 
tem+reservoir is given by H s+R = H(t) + H R + V, where 
H(t) is periodic with period T and the interaction has 
the form V = X €5 Y, where the operator X acts on the 
environment, and Y on the system. 



The Master equation for the reduced density matrix of 
the system can be presented in the forrrPSl 



j kV + ^r^ ifcV )]p 7 ,(t)e < ^.''-^.-*')* 



7<5;&:,/c' 



(ii) 



having defined 

r+ 



with 



1 a0j8,kk' 



^aP,k 



— ^Yap^Y-fS.k'l^S.k' 

= ^5^{<s> a {t)\Y\<s> p {t)y™dt 

= C(X(±t") X) eM-i^,kt")dt" 



^a/3,k 



fa 



As it is evident from Eq. (11), Floquet theory allows 



one to treat the time periodic case with a formalism 
which is formally identical to the one used in the time- 
independent case. The relevant effects due to the periodic 
driving are captured by the use of the Floquet basis. 

Further simplifications can be made if the secular ap- 
proximation holds^as we will assume in the rest of the 
paper). As in the time- independent case, the equations 
for the populations decouple from those for the (off- 
diagonal) coherences. In the steady state the coherences 
vanish. The populations are given by a "detailed bal- 
ance" condition 



St 

act 



(12) 



Ws^c = i E \ Y ^k\ 2 9(&s a ,k) 



and 



+00 



(X(t) X)e iu)t dt . 



(13) 



(14) 



This property of the stationary solution, when expressed 
in the Floquet representation, is crucial to obtain the 
pumped charge also in the dissipative case. 



IV. FLOQUET APPROACH TO COOPER PAIR 
PUMPING 



Equipped with the Floquet formalism outlined above, 
we now derive an expression for the pumped charge both 
in the absence and in the presence of an external envi- 
ronment. We first consider the case of a unitary evolu- 
tion where, in the adiabatic case, we should recover the 
known relation, Eq. ^ , between pumping and geometric 
phases. 
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A. Pumped charge in a closed system 

At first we ignore any coupling with the external en- 
vironment. The dynamics is unitary and governed by a 
time-periodic Hamiltonian. It is meaningful to compute 
the pumped charge in a given cycle only for those states 
that, up to a phase, do come back to their initial value. 
These are the Floquet states (one should keep in mind 
that no assumption of an adiabatic dynamics is done at 
this point). By employing the Schrodinger equation in 
Eq.fJlJ it is possible to write it in the form 



Q 



(tr) _ 



-2ie 



d t (*(t)\d v Mt))dt 



(15) 



Since the average is performed over a Floquet state, de- 
fined in Eq.^, it straightforward to obtain 



Q(.*r) = _ 2e |^e f » . 



(16) 



Eq.(16l is the first result of this paper. It gives a gen- 
eral formula for the transferred charge in a supercon- 
ducting circuit which is valid both under adiabatic and 
non-adiabatic conditions. Obviously it reduces to Eq.Q 
in the adiabatic case. 

The pumped contribution to the transferred charge can 
be obtained by subtracting from Eq.( 16 1 the supercurrent 



term (associated to the dynamical phase) 

Q p = Aire k d v {c a , k (Lp)\c a .k{tp)) 



(17) 



B. Pumped charge of a superconducting circuit 
coupled to an external environment 

The Floquet approach allows for a very simple and ap- 
pealing generalization of Eq.(17) to the dissipative case. 



For simplicity we consider the case in which the secu- 
lar approximation holds. The transferred charge for a 
system defined by a density matrix is given by 



)(*>-) _ 



2e 



Tr 



((V?(i,v?))p st (^))di (18) 



where p st is the reduced density matrix of the system in 
the steady state (we are interested in obtaining a suitable 
expression for the pumped charge in the stationary limit 
after all transient effects have disappeared). Noting that 
(see the previous Section) in the Floquet basis all the 
coherences vanish in the long time limit, the expression 
for the transferred charge takes the form 



(19) 



Hence, the charge passing through the circuit is the 
weighted average of the charge which would have passed 



if the system had been in a pure Floquet state, see 
Eq.(16|. The weights are the populations of these Flo- 
quet states in the quasi-stationary case. The previous 
expression can be split in a geometric and a dynamic 
part and the pumped charge is given by 



4neJ2k P s U<P)d v 

k . v 



(20) 



This is the central result of our work. Eq.(20) reduces 
to all known cases in the corresponding limits. In ad- 
dition allows us to explore regimes that have not been 
considered so far. The form in the dissipative case is 
self-explaining, it is the average of the corresponding ex- 
pression in the noiseless case weighted by the populations 
of the Floquet states. 

In t he r est of the paper we will apply the general result 
of Eq.(20) to the Cooper pair sluice^ which, as we will 



see in the next section, can be described by a two-level 
Hamiltonian. It is therefore useful to give explicit formu- 
las in the case of a two-dimensional Hilbert space. In this 
case we have only two independent Floquet states that 
we call \^ a ) and |^s)- The populations in the stationary 
state are given by 



p st = 

r aa 




w p . 




Wp. 


+0, + 




Pup - 




w a . 


-+(8 


Wp. 


+ a + 





(21) 



Replacing this in Eq.( 20 ) and exploiting the relation e a — 
—tp we obtain 

Qp = 4ire^f — r ^ ■ w — kd v {c a ^{(p)\c a ^KV)) ■ 



(22) 



V. PUMPING IN THE COOPER PAIR SLUICE 

The Cooper pair sluice^ is a superconducting transis- 
tor where the pumping effect is achieved by a modula- 
tion of the Josephson couplings and the gate voltages. A 
sketch of the sluice is shown in Figj2j The central island 
is connected to the two superconducting leads by two 
tunable Josephson junctions. Its charging energy can be 
tuned by means of a gate voltage. The Hamiltonian of 
the sluice is 



H = E c (h 



^2 J i C0S ( 
=L/R 



(23) 



where Ec is the charging energy of the central island, Jj 
are the Josephson coupling to the left (i = L) and right 
(i = R) electrodes and n g is the gate charge which can 
be modulated by changing the gate voltage V g as shown 
in Fig[2] In order to tune the Josephson couplings the 
junctions are replaced by SQUIDs which behave as sin- 
gle Josephson junctions with an effective coupling which 
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FIG. 2: Circuit scheme of the Cooper pair sluice. 





FIG. 4: Pumped charge vs. in the adiabatic case computed 
with the Floquet theory (red crosses) confronted with the 
same quantity computed with the analytical expression (241 
(black solid line), < <p < 2ir. It is Jmin/Jmax = 0.03, J m ax = 
O.lEc, TEc/h = 8400. Charge is in units of e. 



two-level Hamiltonian 
1 



H 



[B x a x + B y a y + B z u z 



(25) 



FIG. 3: The variation in a cycle of n g , Jl e Jr. These latter 
are expressed in units of Ec- In this graph n g varies between 
n g nun = 0.2 and n g max = 0.8, Jl and Jr vary between 



iTrr 



0.003£ c and J n 



Q.1E C . 



can be varied by changing the flux piercing the loop 
J($) = j(°) cos(7T$/$ ) (with $ the flux through the 
loop and $o the flux quantum). The charge h on the 
central island and the phase 9 are canonically conjugated 
variables. The Hamiltonian is varied along a cyclic path, 
we change periodically $l, Qr e V g , determining a peri- 
odic variation of n g , Jl e Jr. In all the cases considered 
here, these parameters are assumed to depend on time in 
the same way as in Pekola et aP^. The time-dependence 
of the parameters is shown in Fig[3j Notice that when J L 
is maximum Jr is minimum and vice versa, meaning that 
when a SQUID is open the other is closed and vice versa. 
In an ideal situation where the minimum value of the 
Josephson couplings could be reduced to zero the charge 
passing through the system in one cycle would be ex- 
actly quantized in units of 2e (the supercurrent contribu- 
tion vanishes in this setup) . In a realistic case where the 
SQUID loops do not close perfectly the pumped charge 
is given by^ 



-2e 1 - 2- 



J„ 



COS (f 



(24) 



If n g during the cycle stays close enough to the degen- 
eracy value 1/2, we can describe the system using the 



where B x = Eq{1 — 2rt g ), B y = J^sini/j, and B z — 
{Jl + Jr cos tp). 

In the following sections we will use the two-level ap- 
proximation to present our results for the pumped charge 
in the Cooper pair sluice. 



1. Unitary evolution 

We first consider the case in which the environment is 
absent. We compute numerically the Floquet exponents 
for the Hamiltonian defined by Eq. ( 25 ) and by means of 
Eq.(22) we obtain the pumped charge. The first case we 



consider is the adiabatic limit to compare our approach 
with the known results. The behavior of the pumped 
charge as a function of the phase bias is shown in Fig|4j 
The results of the Floquet approach are tested against 
the analytic resullP, Eq.@. We chose TE c /h = 8400, 



a value which guarantees amply to be in the adiabatic 
limit. The reason for this large value is due only to 
the simplicity to compute numerically the Floquet ex- 
ponents. Obviously we do not expect any changes in the 
results as long as we are in the adiabatic regime. The 
numerical calculations agree well with the analytical ex- 
pression of Eq. ( 24 ) . 

As we already discussed, the Floquet approach to 
pumping allows us to go beyond the adiabatic regime. 
An example is shown in Fig js] for TEc/h = 2.1. In this 
particular case the pumped charge is much smaller than 
that obtained in the adiabatic regime. This example was 
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J»in/Bc = 
Jmin/Ec =0.003 
Jrain/Ec =0.006 
Jmin/gQ =0-009 



FIG. 5: Pumped charge vs tp with TEc/h — 2.1 for different 
values of J m i n , < tp < 2ty. Charge is in units of 10 -5 e. 



indeed chosen just to demonstrate the power of the Flo- 
quet approach. There are cases, for a suitable choice of 
the parameters' loop, in which it is possible to obtain 
charge quantization also under non-adiabatic conditions. 

It is interesting to note that in the non-adiabatic case 
the pumped charge is phase-dependent also in the case in 
which J m in = 0. This case is very similar to the Cooper 
pair shuttlePEZl. 
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FIG. 6: Population p B aa of the Floquet state with the lower 
quasi-energy vs temperature O (in its dimensionless form) for 
different values of tp. We are in the adiabatic limit (TEc/h = 
8400 and Ai5 m ; n ~ Eq/W). Temperature is expressed in 
units of Ec/ks- (Inset) The same graph magnified to show 
the dependence on tp. 



2. Influence of the external environment 

The second application of our pumping formula deals 
with the case in which the sluice is coupled to an external 
environment. The main source of decoherence is due to 
charge fluctua tions. In terms of the coupling Hamitonian 
the operator of the system Y cx a z 
where 5n a = C g V g /2e expresses the 



introduced in III C 
and X = EcSrig, wuao u,u g 
fluctuations of the gate voltage (C g is the gate capaci- 
tance) . As usual, we assume that the charge fluctuations 
are due to the thermal noise of a resistance R! put in 
series with C g . There are also some fluctuations in the 
fluxes $l and <$>r, but these ones are coupled to the 
Josephson energies Jl and Jr which are Jl,Jr "C Eq , 
hence these fluctuations are much smaller than the charge 
fluctuations which are coupled to Ec and we can neglect 
them. The function g{u) = J™ {X(t") X)e luJt " At" which 



appears in Eq.(13l is defined as 



( coth 



(*r) 



(26) 



In all the subsequent calculations R = 300kSl. 

There are several methods to compute the Floquet 
quasi-spectrurrP^l. We diagonalized numerically the Flo- 
quet operator H{t) in the composite Hilbert space^ as 
briefly discussed in Section |III A| From the knowledge 
of the Floquet eigenvalues and eigenvectors it is possible 



to compute the steady state populations in each Floquet 
mode. These latter are shown in Figj6j The parameters 
are chosen to be such that the pumping is adiabatic. 



By means of Eq.(22| the pumped charge in a dissipa- 
tive case is readily obtained. An example is shown in 

As it should be expected, on increasing the tempera- 
ture the pumped charge decreases. The pumped charge 
associated to the two Floquet states is opposite for a 
given cycle. Therefore the progressive mixture of the 
two states suppresses the size of the pumping. 



3. Floquet states in the adiabatic limit 

In the limiting case of an adiabatic cycle our approach 
to pumping should reduce to the one discussed in Ref . [27] 
In the remaining of this section we will address this point 
by analyzing the Floquet states in the adiabatic approx- 
imation. Floquet states are the eigenstates of the op- 
erator H(t) — ihdt- We assume a non-degenerate spec- 
trum with instantaneous eigenvectors and eigenvalues 
given respectively by {\k(t))} and {Ek(t}}. Adiabatic 
condition require that a = 1/(T min u)ki) <C 1 where 
<jJki{t) = (Ek(t) — Ei(i))/h. In the basis of the instanta- 
neous eigenstates the operator to be diagonalized reads 
E k (t)5 kl - hw kl (t) where w kt (t) = i(k(t)\i(t)) . To first 
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ip 


= 0.8 


If 


= 1.6 


ip 


= 2.3 
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FIG. 7: Modulus of the pumped charge vs temperature for 
different values of tp. We are in the adiabatic limit (TEc /h = 
8400 and A£ m j n — Ec/W). Temperature is expressed in 
units of Ec/kB and charge in units of e. 



order in a the corresponding eigenvectors are 



\$ k{t )) = \k(t))-^^\l(t)) + 0(a 2 ) 



l^k 



(27) 



with eigenvalues 

cT 



^ (£ fc (t)-/«« fcfe (i))di + 0(a 2 ). 



(28) 



It can be shown that these eigenvalues are invariant under 
gauge transformations. These states are eigenstates of 
the operator 



where D is the the transformation from a given (time- 
independent) basis of the system Hamiltonian. It was 
argued in Ref. [57] that the system relaxes in the basis of 
Eq. |29| ) . This is what is contained in the result of Eq. ( 17 ) 
for the pumped charge. 



VI. CONCLUSIONS 

In this paper we developed a new approach to Cooper 
pair pumping using Floquet theory of periodically driven 
quantum systems. We found that the pumped charge 
can be expressed in a very natural way in terms of the 
Floquet eigenstates (in the adiabatic limit there is a 
clear connection between Floquet exponents and states 
and the geometric phase acquired by the system). This 
approach does not require the adiabatic limit to hold, 
and can be used to work out the pumped charge out of 
the adiabatic regime (see Fig (5]). We further extended to 
the dissipative case where we provided a general formula 
to compute pumping. In order to demonstrate the power 
of our approach we applied it to the case of Cooper 
pair sluice. In the known limits we recovered previous 
results. We further discussed new regimes which are 
now addressable due to the expressions for the pumped 
charge given in Eq. ( 20 1 . 



D(tyH(t)D{t) - ihD(tyD{t) 



(29) 



Acknowledgments 

We acknowledge very useful discussions with A. Shnir- 
mann and J. P. Pekola. This work was supported by EU 
through the projects GEOMDISS, QNEMS, SOLID and 
NANOCTM. 



1 D.J. Thouless, Phys. Rev. B 27, 6083 (1983). 

2 P.W. Brouwer, Phys. Rev. B 58, R10135 (1998). 

3 L.P. Kouwenhoven, A.T. Johnson, N.C. van der Vaart, 
C.J.P.M Harmans, and C.T. Foxon, Phys. Rev. Lett. 67, 
1626 (1991). 

4 H. Pothier, P. Lafarge, C. Urbina, D. Esteve, and M. H. 
Devoret, Europhys. Lett. 17, 249 (1992). 

5 F. Zhou, B. Spivak, and B. Altshuler, Phys. Rev. Lett. 82, 
608 (1999); Yu. Makhlin and A. D. Mirlin, ibid. 87, 276803 
(2001); O. Entin-Wohlman, A. Aharony, and Y. Levinson, 
Phys. Rev. B 65, 195411 (2002); M. Moskalets and M. 
Biittiker, Phys. Rev. B 66, 205320 (2002); I. L. Aleiner 
and A. V. Andreev, Phys. Rev. Lett. 81, 1286 (1998); M. 
Blaauboer and E. J. Heller, Phys. Rev. B 64, 241301 (R) 
(2001); B. L. Hazelzet, M. R.Wegewijs, T. H. Stoof, and 
Yu. V. Nazarov, ibid. 63, 165313 (2001); R. Citro, N. An- 
drei, and Q. Niu, ibid. 68, 165312 (2003); P. W. Brouwer, 
A. Lamacraft, and K. Flensberg, ibid. 72, 075316 (2005); 



L. Arrachea, A. Levy Yeyati, and A. Martin-Rodero, ibid. 
77, 165326 (2008); A. R. Hernandez, F. A. Pinheiro, C. 
H. Lewenkopf, E. R. Mucciolo, Phys. Rev. B 80, 115311 
(2009); J. Splettstoesser, M. Governale, J. Konig, and R. 
Fazio, Phys. Rev. Lett. 95, 246803 (2005); E. Sela and Y. 
Oreg, ibid. 96 166802 (2006); D. Fioretto and A. Silva, 
ibid. 100, 236803 (2008). 

6 L. J. Geerligs, S. M. Verbrugh, P. Hadley, J. E. Mooij, 
H. Pothier, P. Lafarge, C. Urbina, D. Esteve, and M. H. 
Devoret, Z. Phys. B: Condens. Matter 85, 349 (1991). 

7 J. P. Pekola, J. J. Toppari, M. Aunola, M. T. Savolainen, 
and D. V. Averin, Phys. Rev. B 60, 9931(1999). 

8 J. P. Pekola and J. J. Toppari, Phys. Rev. B 64, 172509 
(2001). 

9 M. Aunola and J. J. Toppari, Phys. Rev. B 68, 020502(R) 
(2003). 

R. Fazio, F.W.J. Hekking, and J.P Pekola, Phys. Rev. B 
68, 054510 (2003). 



9 



11 A. O. Niskanen, J. P. Pekola, and H. Seppa, Phys. Rev. 
Lett. 91, 177003 (2003) 

12 M. Governale, F. Taddei, R. Fazio and F. W. J. Hekking, 
Phys. Rev. Lett. 95, 256801 (2005). 

13 M. Mottonen, J. P. Pekola, J.J. Vartiainen, V. Brosco, and 
F.W. J. Hekking, Phys. Rev. B 73, 214523 (2006). 

14 R. Leone, L. P. Levy, and P. Lafarge, Phys. Rev. Lett. 100, 
117001 (2008). 

15 V. Brosco, R.Fazio, F.W.J. Hekking and A. Joye, Phys. 
Rev. Lett. 100, 027002 (2008). 

16 J.-M. Pirkkalainen, P. Solinas, J. P. Pekola, and M. 
Mottonen, Phys. Rev. B 81, 174506 (2010). 

17 J. J. Vartiainen, M. Mottonen, and J. P. Pekola, Appl. 
Phys. Lett. 90, 082102 (2007). 

18 M. Mottonen, J.J. Vartiainen, and J. P. Pekola, Phys. Rev. 
Lett. 100, 177201 (2008). 

19 J. E. Avron, A. Elgart, G.M. Graf, L. Sadun, Phys. Rev. 
B 62, R10618 (2000). 

20 H.-Q. Zhou, S. Y. Cho, and R. H. McKenzie, Phys. Rev. 
Lett. 91, 186803 (2003). 

21 G. Falci, R. Fazio, G. M. Palma, J. Siewert, V. Vedral, 
Nature 407, 355 (2000). 

22 W. Xiang-bin, M. Keiji, Phys. Rev. B 65, 172508 (2002). 

23 A. Blais, A. M. S. Tremblay, Phys. Rev. A 67, 012308 
(2003). 

24 P. J. Leek, J. M. Fink, A. Blais, R. Bianchetti, M. Gppl, J. 
M. Gambetta, D. I. Schuster, L. Frunzio, R. J. Schoelkopf, 
A. Wallraff, Science 318, 1889 (2007). 

25 J. Jones, V. Vedral, A. Ekert, and G. Castagnoli, Nature 
(London) 403, 869 (2000). 

26 P. Zanardi, M. Rasetti, Phys. Lett. A 264, 94 (1999). 

27 J. P. Pekola, V. Brosco, M. Mottonen, P. Solinas, and A. 



Shnirman, Phys. Rev. Lett. 105, 030401 (2010). 

28 P. Solinas, M. M ottonen, J. Salmilehto, and J. P. Pekola, 
larXiv: 1007.23471 

29 J. Salmilehto, P. Solinas, J. Ankerhold, and M. Mottonen, 
Phys. Rev. A 82, 062112 (2010). 

30 E. Sjoqvist, A. K. Pati, A. Ekert, J. S. Anandan, M. Er- 
icsson, D. K. L. Oi, and V. Vedral, Phys. Rev. Lett. 85, 
2845 (2000); I. Fuentes-Guridi, J. Pachos, S. Bose, V. Ve- 
dral, and S. Choi, Phys. Rev. A 66, 022102 (2002); G. 
De Chiara and G. M. Palma, Phys. Rev. Lett. 91, 090404 
(2003); A. Carollo, I. Fuentes-Guridi, M. F. Santos, and V. 
Vedral, ibid. 90, 160402 (2003); M. S. Sarandy and D. A. 
Lidar, ibid. 95, 250503 (2005); R. S. Whitney, Y. Makhlin, 
A. Shnirman, and Y. Gefen, 94, 070407 (2005); P. Thun- 
strom, J. Aberg, and E. Sjoqvist, Phys. Rev. A 72, 022328 
(2005); G. Florio, P. Facchi, R. Fazio, V. Giovannetti, and 
S. Pascazio, Phys. Rev. A 73, 022327 (2006); D. Parodi, 
M. Sassetti, P. Solinas, P. Zanardi, and N. Zanghi, ibid. 
73, 052304 (2006). 

31 M. Wubs, K. Saito, S. Kohler, P. H&nggi, and Y. 
Kayanuma, Phys. Rev. Lett. 97, 200404 (2006). 

32 M. Moskalets and M. Biittiker, Phys. Rev. B 66, 205320 

(2002) ; ibid. 70, 245305 (2004). 

33 M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 

34 S. Guerin and H. R. Jauslin, Adv. Chem. Phys. 125, 147 

(2003) . 

35 H. Sambe, Phys. Rev. A 7, 2203 (1973). 

36 L. Y. Gorelik, A. Isacsson, Y. M. Galperin, R. I. Shekhter, 
and M. Jonson, Nature 411, 454 (2001). 

37 A. Romito, F. Plastina, and R. Fazio Phys. Rev. B 68, 
140502(R) (2003). 



